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In this paper, we study efficient algorithms towards the construction of any arbitrary Dicke state. 
Our contribution is to use proper symmetric Boolean functions that involve manipulations with 
Krawtchouk polynomials. Deutsch-Jozsa algorithm, Grover algorithm and the parity measurement 
technique are stitched together to devise the complete algorithm. Further, motivated by the work 
of Childs et al (2002), we explore how one can plug the biased Hadamard transformation in our 
strategy. Our work improves the results of Childs et al (2002). 
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I. INTRODUCTION 

Multipartite entanglement is one of the important ar- 
eas in the field of quantum information that has many 
applications including quantum secret sharing. In this 
paper, we focus on the Dicke states which are useful 
building blocks in realizing multipartite entanglement. 
The n-qubit weight w Dicke state, is the equal su- 

perposition of all n-q ubit states of weight w. We refer 
to (H, i, d EL El El and the references therein for 
detailed discussion. 

While the main focus from the viewpoint of experimen- 
tal physics is to actually provide the implementation of 
specific Dicke states, our focus is from theoretical algo- 
rithmic angle and the only result presented in this direc- 
tion appeared in In this work, we show how one can 
efficiently construct Dicke states by using the combinato- 
rial properties of symmetric Boolean functions, two well- 
known quantum algorithms, and the generalized parity 
measurement. By efficient, we mean that the resource 
requirements in terms of quantum circuits and number 
of execution steps is poly(n) to obtain \D™). 

Let us consider n-qubit states in the computa- 
tional basis {0, l} ra that can be written in the form 

£xe{o,i}™ a x\%), where X^ G {o,i}" Kl 2 = 1 - Thus ' x can 
also be interpreted as a binary string and the number of 
l's in the string is called the (Hamming) weight of x and 
denoted as wt{x). Based on this an arbitrary Dicke state 
can be expressed as follows: 



Let us also define a symmetric n-qubit state as 



\S n ) = ^2 a wt{x) \x), where 

x£{0,l} n i=0 
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First, we show how one can prepare a symmetric n- 
qubit state with the property that ("^a^l 2 is ^(-^0 
by using Deutsch-Jozsa algorithm Jj|. This requires cer- 
tain novel combinatorial observations related to symmet- 
ric Boolean functions. Then the quantum state out of 
Deutsch-Jozsa algorithm is measured using the parity 
measurement technique Q to obtain |-D") with a proba- 
bility f2(-^). Thus, 0(y/n) runs are sufficient to obtain 
the required Dicke state. Note that a direct approach 
to construct a symmetric state has been presented in [l[ 
using biased Hadamard transform. While the order of 
probability to obtain Dicke state by ours and that of [l[ 
are the same, enumeration results show that the exact 
probability values arc better in our case than that of [l[ . 

Further, motivated by the idea in we improve our 
algorithm further with a modified Deutsch-Jozsa opera- 
tor that involves the biased Hadamard transform. Since 
biased Hadamard transform also helps to generate the 
target symmetric state, the overall probability to obtain 
the Dicke state increases. 

Finally, we can also apply the Grover operator Q be- 
fore the measurement. Since Grover algorithm amplifies 
the amplitude of target symmetric state, this helps to 
reduce the necessary number of steps into O(tfn). 



II. PROPERTIES OF SYMMETRIC BOOLEAN 
FUNCTIONS 

A. Walsh Spectrum of Symmetric Boolean 
Functions 

A Boolean function on n variables may be viewed as 
a mapping from {0,1}" into {0,1}. Let us denote the 
addition operator over GF(2) by ©. Let x = (xi, . . . , x n ) 
and oj = (uji, . . . ,u! n ) both belong to {0,1}™ and the 
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inner product 



x ■ w = x\uji 



' X n Ur, 



Let f(x) be a Boolean function on n variables. Then the 
Walsh transform of f(x) is a real valued function over 
{0, 1}™ which is defined as 



W f (u) 



E 

x€{0,l}™ 



(-!)/(* 



)®x-cj 



An n- variable Boolean function / is called Symmetric 
if f(x) = f(y) for all x, y e {0, 1}" such that wt(x) = 
wt(y). Henceforth, we will denote the set of n- variable 
symmetric Boolean functions as SB n . 

In the truth table of / 6 SB n , it is enough to provide 
outputs corresponding to different weights of elements of 
{0, 1}™ only. So an n- variable symmetric function can be 
expressed by an (n + 1) length bit string as 

re f = L/o,/i, •■•,/«], 

where /j is the output at the inputs of weight i and re / is 
referred to as the simplified value vector. When / g SB n , 
one may note that Wj(x) = Wf(y) for all x, y £ {0, 1}™ 
such that wt(x) = wt(y). Therefore, the Walsh spectrum 
of / can be represented by an (n + l) length integer string 

rwf = [rw f (0) , rw f (1) , . . . , rwf(n)}, 

where rwf(i) represents the Walsh spectrum value at the 
inputs of weight i. 



B. Relation between Walsh Spectrum and 
Krawtchouk polynomials 



We now relate the Walsh spectrum of the symmet- 
ric functions [13] with Krawtchouk polynomials [ill (HI ■ 
Krawtchouk polynomial of degree i is given by 



^( 7? ,n) = ^(-l) J 

3=0 



n — n 



From [13], we get that if wt(ui) = k, then 

71 



i=0 

The (n + 1) x (n + 1) matrix which has Ki(k, n) as the 
i, k)-th element is known as the Krawtchouk matrix [5], 



For example, let us present the Krawtchouk matrix for 
n — 5 and 6 as follows: 
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6 4 2 -2 -4 -6 

15 5 -1 -3 -1 5 15 

20 -4 4 -20 

15 -5 -1 3 -1 -5 15 

6 -4 2 -2 4 -6 

1-11-11-1 1 



In these two matrices, one can verify the properties re- 
lated to the Krawtchouk matrix given in Proposition [TJ 
To determine all the Walsh spectrum values of / € 
SB n , it is enough to multiply ((— 1)^°, . . . , (— 1)^") with 
the (n + 1) x (n+1) Krawtchouk matrix. Applying 
Krawtchouk matrix, the analysis of the Walsh spectra 
of symmetric functions becomes combinatorially inter- 
esting. Elements of a Krawtchouk matrix have nice com- 
binatorial properties and they follow nice symmetry [Toj 
too. We list some of them in the following proposition. 

Proposition 1 1. Ko(k, n) = 1, Ki(k, n) = n — 2k, 

2. (i + l)K i+1 (k,n) = (n-2k)Ki(k,n) 
-{n - i + l)A' l _i ! „(fc, n), 

3. Ki(k,n) = (-l) k K n _i(k,n), 
I ($K i {k,n) = (»)K k (i,n), 

5. Ki(k,n) = (-l)*Ki{n -k,n), 

6. (n-k)Ki(k + l,n) = (n-2i)K z {k,n)-kK t (k-l,n), 

7. (n-i + l)Ki(k,n + l) = (3n - 2i - 2k + l)A,(fc, n) 
-2(n- k)Ki{k,n- 1). 



C. Implementation of Symmetric Boolean 
Functions 

The symmetric Boolean functions can be efficiently im- 
plemented. As described in Q, the circuit complexity of 
n- variable symmetric Boolean functions is 4.5n + o(n). It 
is known that given a classical circuit /, there is a quan- 
tum circuit of comparable efficiency which computes the 
transformation Uf that takes input like \x,y) and pro- 
duces output like \x,y © f(x)}. Thus, we will consider 
that for / G SB n , the quantum circuit U / can be effi- 
ciently implemented using 0{n) circuit complexity. 



III. ALGORITHMS 

A. Find a Special Symmetric Boolean Function 
which maximizes the Walsh Spectrum 

Consider that we want to maximize the Walsh spec- 
trum value corresponding to weight w points and nat- 
urally, from the property of symmetric functions, all of 
them will be equal. Now we present an important com- 
binatorial result to show how to find such symmetric 
Boolean functions. 

Theorem 1 Consider f G SB n . The function f , repre- 
sented as rej, for which the Walsh spectrum correspond- 
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ing to the w weight points will be maximized, can be writ- 
ten as 

!0, ifKi(w,n)>0 
1, ifKi{w,n)<0 
Oorl, ifKi(w,n)=0 

Proof: We have W f {uj) = ELot- 1 )^^^")- 
One may note that the maximum value of 

EHoC- 1 )-^^) is J2to\ K ^ n )\- This is at- 
tained when we take the function of the form as 
described in the theorem. ■ 

Example 1 As example, consider n = 6,w = k = 2. In 

the corresponding column of the (6 + 1) x (6 + 1) matrix, 
we get the values as 1,2,-1,-4,-1,2,1. Thus, we will 
consider the function with ref as [0,0,1,1,1,0,0]. For 
such an f G SBq, the Walsh spectrum values at the points 
uj, such that w = wt(cu) = 2, will be maximized, which is 
1 + 2 + 1 + 4+1 + 2 + 1 = 12. 



putting i — 21. Hence, 



i=0 



Now let us prove this for odd n 

i— 

2 



For n odd and k = ^ip-, from Proposition [T]^2) we get 



,n-l 1 ,n- 1 

A,-+i( ,n) = A, , n) 

+n 2 ' i+l v 2 ' 

n — i + 1 .71—1 . 

Ai_i( , n). 

i+1 v 2 ' ; 

That is, the recurrence relation is as follows: 

t, , n ~ 1 \ 1„ ,71 — 1 , 

^(-j-.n) = - A *-i(^-'") 

77 — 7 + 2 r ^ .77,-1 

: Ki-2{— 

% 2 

One can now show by induction that 



B. Walsh Spectrum of the Special Symmetric 
Boolean function by Combinatorial Property of 
Krawtchouk Matrix 

Next we present certain results related to column sum 
of Krawtchouk matrix. 

Lemma 1 E^o I = £?=o l*i(Li J , «)l = 



,71 — 1 . . ,71 — 1 n — 1 

A 2 i(— ^— ,nj = A 2 i+i(— ^— Vi, 1 <i< 

Using the above two identities and induction, one can 
verify that \K 2 i(^,n)\ = ("y). Thus ; 



j2\m^,n)\=2Y; 



2 /n-1 
2 



Proof: Let us first prove this for even n. 
Following Proposition QJ2), we have 

(i + l)K l+1 (k,n) = (n-2k)Ki(k,n) 

— (n - i + l)A~j_ lin (fc,n). 

For n even, and k = %, we get, 



where i = 21. Hence, we get, 



^|A l (L^j,n)|=2.2^ i =2- 



Using Proposition [TJ5) , we get that 



Ki+x{-,n) = — — Ai_i )n (-,n). 

2 l + 1 2 

That is, the recurrence relation follows: 



n n-i+2 n 

^,(-,77) = : A i _ 2 ,„(-,n), 

2z2 



with the initial conditions A"o(^,n) = 1 and A 1 (^,n) = [^] or [_§ J ■ T7i 
as available from Proposition QJl). Thus one may note 
that for odd i, Aj(^,n) = 0. Further, using induction, 
for even i, we get 



^|A,(L-J,n)H^|A l (r-l,n)|. 

That completes the proof. ■ 

Theorem 2 Let f G <S£>„ 6e as explained in Theorem[J\ 

towards maximizing the Walsh spectrum values at weight 



;i)(~'<rfi'» 2 = ( L fj) ( ™' < ^ ))2is!J< S ) 



Thus, 



i=0 



2 1/ " \L2- 

Proof: The Walsh spectrum in this case is 



rw f (\-]) = ^|A-(r-],7,)| 

2=0 

n 

= ^|A-(^j,77)|. 

i=0 



4 



c(n) 



FIG. 1: Plot of c(n) vs n for 1 < n < 100. 



Thus the total sum of the squares of the Walsh spectrum 
values at weight \^~\ or [f J is 



rti 



(Ei^(^i^)i) 



which is f2(2-;=), by Stirling's approximation. ■ 
One may similarly note that for the trivial cases of 
w = or n, if one chooses / £ SB n following Theorem [TJ 
then (™) (rwf(w)) 2 = 2 2n . However, proving the result 
similar to Theorem [5] for any n and any weight w, in gen- 
eral, seems to be quite tedious. Thus we make detailed 

• / \ • n ( " )( rw f( w )) 2 

enumerations to obtain c(n) = mm^^Q — that 

has been verified for n < 1000 and we note that the val- 
ues stabilize as c(999) = 1.24793 and c(1000) = 0.797685. 
The graph of this is plotted in Figure [TJ for 1 < n < 100, 
the points for odd n are coming above and those for even 
n are coming below. Since we are not providing a proof 
of this, we refer this as follows. 

Fact 1 Let f £ SB„ be as explained in Theorem Q] to- 
wards maximizing the Walsh spectrum values at weight 
w. Then the total sum of the squares of the Walsh spec- 
trum values at weight w, (™) (rwf(w)) , is fl^ 2 -^). 



C. Relation between Deutsch-Jozsa algorithm and 
the Walsh Spectrum of Symmetric Boolean Function 

Given / is either constant or balanced, if the cor- 
responding quantum implementation Uf is available, 
Deutsch-Jozsa Q provided a quantum algorithm that de- 
cide in constant time which one it is. Let us now describe 
our interpretation of Deutsch-Jozsa algorithm in terms 
of Walsh spectrum values. We denote the operator for 
Deutsch-Jozsa algorithm as 



where the Boolean function / is available as an oracle U / . 
For brevity, we abuse the notation and do not write the 
auxiliary qubit, i.e., ^^j 1 ^ and the corresponding output 
in this case. 



Now one can observe that UA 



»/|0> 



'Sin 



E E 

ze{o,i}" xt{o,i}" 



(-1)*- 



B/(z) 



\z) 



zG{0,l}" 



2" 



Note that the associated probability with a state \z) is 

W ^}^ . Hence we have the following technical result as 
pointed out in 14| with our interpretation for symmetric 
functions. 

Proposition 2 Given an n-variable Boolean function f , 
P/|0)®" produces a superposition of all states z £ {0, 1}™ 

with the amplitude corresponding to each state \z) . 

Specially, if f £ SB n , then the amplitude corresponding 
to \z) z S SLaCzlteM^) , 



D. Algorithm 1: Deutsch-Jozsa Algorithm with 
Special Symmetric Function 

Based on the overall properties, we propose a quantum 
algorithm as shown in the following algorithm. 

Algorithm 1 

1. Choose / £ SB n as explained in Theorem Q] to max- 
imize the Walsh spectrum values at weight w. 

2. Use the Deutsch-Jozsa algorithm to obtain a sym- 
metric n-qubit state \S n ) = a wt(x)\ x )> such that 

xe{o,i}" 

OKI 2 iflfi(^). 

3. Apply the parity measurement strategy. If the an- 
cilla state is measured at the basis U w \(), then |-D") is 
successfully obtained. Else go to Step 2 and iterate. 

The following result provides the estimate of complex- 
ity of our algorithm. 

Theorem 3 Let f £ SB n be as explained in Theorem Q] 
towards maximizing the Walsh spectrum values at weight 
w. Given Fact [3 Deutsch-Jozsa algorithm produces a 
symmetric n-qubit state (before the measurement) \S n ) = 
E a w t(x)\x), such that (™)|aj 2 is ^(7^)- 



ce{o,i}' 



V, 



H® n U f H 



Proof: The proof follows from Theorem [TJ Fact [TJ and 
Proposition O ■ 
Now the final step is to measure the symmetric state 
until we get the target Dicke state by using parity mea- 
surement method [8j, Section IIIA] . Now we explain how 
to exploit the parity measurement in our purpose. Note 
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FIG. 2: Generalized Parity Module as in g, Fig. 1]. If 
wt(x) = w, then the ancilla will become U W \Q. 



P, Section IIIA] assumes to use n dimensional qudit an- 
cilla, but we consider a qudit |£) of dimension n + 1 here. 
A certain unitary operator U is designed such that, 

ICc/ICf/'IC.-.^'^IC^IC) 

are all orthogonal to each other and U n+1 \(,) = Since 
|C) is an (n+ l)-dimensional state, one can indeed obtain 
a set of such n + 1 orthogonal states. The parity mea- 
surement is done on the 

{|C),c/|C),t/ 2 IC),...,f/"IC)} 

basis. Here \Q is used as the target state. For the n-qubit 
control state \x), if it has weight w then its corresponding 
target state, after application of this circuit, will become 
U W \Q (see Figure [2]). Now consider an n-qubit symmetric 
state as the control input, which is |r) = a w \x), 

xe{o,i}" 

where w = wt(x), J27=o (i)l a *| 2 = 1- After applying this 
circuit, one obtains a w \x)U w \0. Thus, if one 

a~e{o,i}™ 

measures in 

{|C>,c/|C>,c/ 2 IC),...,f/™IC>} 

basis, then the state |-D") will be obtained when the mea- 
surement result of the ancilla state is U W \Q. 

Since the probability of target Dicke state is f2(-^=), 

we should repeat the whole procedure at most 0{y/n). 

Example 2 Let us have an example taking n = 6, w = 2 
to outline our method. In this case the Dicke state will 
be 



\D* 



E 

x£{0,l} 6 ,wt{x)=2 



1 



We start with is an n — 6 variable symmetric Boolean 
function having Walsh spectrum value at each of the 
weight w = 2 point as 12 (following Theorem [71 one 
may refer to Example{l\ also). There are (™) = (2) = 
15 such points. The amplitude associated with each 

we 



= -i . Thus 

16 



point after the DJ algorithm is |f 

9 e t Eig{0,l} e ,K.((i)=2 T&\ X ) +Si€{0,l} 6 ,iot(i^2 fl il 3: ) ml ~ 

tially. Thus, the probability associated with \D^} will be 
GKtI) 2 = Ml = 0.52734375 and hence one may note 
that the probability that, after parity measurement, it will 
land into \D®) is quite high. 



Prob- 
ability 




FIG. 3: Plot of probabilities associated with \D^) against w 
in our case (above) and in [l|] (below) for n — 999. 



E. Comparison with a Previous Method 



It was explained in |8| how one can obtain | D^) from 
-\r X^efo 1}™ \ x )- However, the idea explained in |8|, Sec- 
tion IIIA] works efficiently only for w = [§] or [-§ J • The 
most general work in this direction has appeared in [l| 
where biased Hadamard transformation was exploited. 
The strategy of [J uses biased Hadamard transformation 



W M W. J 

n V nj 

on \Q)® n such that the probability associated with \D™) 

will be ( n )(H)^(i _ w\n-w > rr i e Qf 1 \ Thus 

the probability of our case and also in [l| are of the same 
order. While the theoretical comparison of the exact 
probability values seems elusive, we have made detailed 
enumerations to observe that the exact probability values 
in our case are better than that of [l| . Below we present 
enumeration results towards that. 

First we present two graphs to show the probability 
values associated to |-D") for all w, when n = 999 (to 
represent odd case) or 1000 (to represent even case). For 

our case, it is ( n ) ( -^-- i = ^ A l } w ' n ^ (after application of 
Deutsch-Jozsa algorithm without measurement), and for 
the case of it is (^)(^) w (l - f ) n ~ w (after applica- 
tion of biased Hadamard transform) . From Figure [3] and 
Figure 21 it is clear to note that our method provides 
higher probability (the upper curve) in all the cases ex- 
cept w = 0, n (which are trivial ones) and w = ^ for 
n = 1000. 

In both figures, the present algorithm shows some vari- 
ation of probability when the weight is around . To 
check whether or not these cases still shows the higher 
probability than the previous method, we look into the 
w = [j\ case a little bit more. From Figure 03 one may 



6 



Prab- 




FIG. 4: Plot of probabilities associated with \D™) against w 
in our case (above) and in 11] (below) for n = 1000. 
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FIG. 5: Plot of probabilities associated with \D^) against n 
in our case (above) and in [l| (below) for n = 4 to 1000 and 

™=LtJ- 



note that our probabilityvalues (the upper curve) are 
higher than the case of pj. These results explains the 
advantage of the use of a suitable symmetric Boolean 
function which shows higher Walsh spectrum values for 
the given weight. 



F. Algorithm 2: Additional Improvement by 
exploiting biased Hadamard Operator 

We have provided numerical evidences that using 
proper symmetric Boolean functions in Deutsch- Jozsa al- 
gorithm provides better probability than the use of biased 
Hadamard transform as described in However, mo- 
tivated by [l|, a natural extension should be to couple 
biased Hadamard transform in Deutsch- Jozsa algorithm 



instead of (unbiased) Hadamard transform. Thus, let 
us refer to the general description of a Hadamard type 
transformation (biased or unbiased) that can be written 



as B Tn 



We will replace the 



standard notation of w here by r as we will not restrict 
ourselves to integer values w <G [0,...,n], but use any 
real number r £E [0, n] to obtain the optimum probability 
of success to get a Dicke state. 

Instead of using the operator V } = H® n U f H® n , let 
us first describe the most general operator of the form 



V' s = Bf^UfB, 



V2 ,m 



(i) 



where r\ 1 r2 are real numbers in [0,n]. 

First, we consider the case when r\ = ^, i.e., B ritn 
H, but T2 = r varies towards optimization. That is 



V' f = H® n U f B* 



(2) 



One may note that the application of V f on |0)® n will 
produce 

^/|0) 8 " = ^E, e{0l l}n 

where d(x, z) is the (Hamming) distance between two 
same length binary strings x and z. Before proceeding 
further, we also have the following technical result. 

Proposition 3 Let V' f = H® n U f Bf^. Iff e SB„ then 
X>'y>|0)®" is a symmetric state. 



f 

Proof: We need 



to prove that 
2 (-) 2 is same 
for all the z £ {0, 1}™ having the same Hamming 
weight. Let us consider u,v <E {0,1}™ such that 
u 7^ v, but wt(u) = wt(v). Then we need to prove 
that £* 6{ o,i}4-ir" ffi/(a:) (l - r)^^(r)^ = 

The proof follows from the fact that 



(-) 2 = 
2 (-L) 2 



V f—-]\x-v®f(x)(-i _ r\2 

Z-ixe{0,l}",wt(x)=w\ L ) \ L n > 

given that / is symmetric. ■ 
The main problem in this case is that we need to go 
for trial and error by modifying the symmetric Boolean 
functions and trying out different values of So 
far, we could not obtain the exact characterization of 
symmetric functions while biased Hadamard transform 
is used. Based on this, we propose Algorithm 2 as follows. 

Algorithm 2 

1. Apply V f to |0)® ra to obtain a symmetric n-qubit 
state. The value of £ and the choice of the symmetric 
Boolean functions are achieved heuristically. 

2. Use parity measurement strategy. If the ancilla state 
is measured at the basis U W \Q, then \D™) is successfully 
obtained. Else take the parameters as in Step 1 once 
again and iterate Step 2. 
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As we could not characterize this, to provide some ex- 
perimental results in this direction (see Table [IJ, we used 
the following method for some small values of n (n = 4 
to 9). We select each of the Boolean functions / from 
SB n . Given / and a specific weight w, 1 < w < n, we 
write the expression of success probability as a function 
of r. Then we apply the Maximize function available 
in Mathematica 8.0 to compute the optimum value of r 
given /, w, so that the success probability becomes maxi- 
mum. Note that, as we could not characterize the process 
yet, this is an exhaustive task and for each n, it requires 
checking of 2 n+1 symmetric Boolean functions. That is 
the reason, we can study it for only a few small values of 
n. However, this is a classical computation that can be 
done as an off-line work. Once such programs are exe- 
cuted, we can have a database of proper / G SB n and the 
corresponding r to have the optimal success probability 
to obtain \D™). Given these data, the actual quantum 
algorithm to obtain Dicke states can be efficiently imple- 
mented. 



IV. NUMERICAL COMPARISON OF THREE 
APPROACHES 

Now we compare three approaches: 

• using biased Hadamard operator as shown in 

• Algorithm 1 based on Z>/|0)®", and 

• Algorithm 2 based on V f\Q)® n . 

The first two cases need 0(y/n) complexity, and the third 
one is a heuristic that shows improved results than the 
first two. Some numerical results of the probability asso- 
ciated with l-D") are shown in Table J for n = 4, ■ ■ • , 9. 
As shown in the table, we note that Algorithm 2 provides 
the highest probability than others. There are a few in- 
teresting observations from the enumeration results. 

• We note that the improvements using V f are 
highly significant at w = |_§J or an d the signif- 
icance reduces as w moves away from the middle, 
i.e., towards w = 1 or n — 1. 

• In case of using 2?'/, the success probabilities at w 
and n — w weights are same for all the values of w, 
i.e., w > 1. However, the values of r in those cases 
are same at w and n — w weights for w > 2 only. 



V. ADDITIONAL IMPROVEMENT BY 
APPLYING GROVER ALGORITHM 

A. Additional Use of Grover Algorithm 

Although we can construct the target Dicke state by 
measuring the intermediate quantum state, we may in- 



crease the efficiency further by using the amplitude am- 
plification method. Based on this, an adiabatic evolution 
has been used towards amplitude amplification of the de- 
sired states in [l[ , but no complexity analysis was shown. 
In this work, instead, we apply the conventional Grover 
algorithm [3] as it provides a quadratic speed-up. 

Instead of equal superposition = H® n \0)® n = 
-4- X)irg{o i}" \ x ) m Grover algorithm, we will use the 
symmetric state of the form |*) = PfflO)®") = 
J2xt{o i}" Wf 2^ I 3 ') exploiting the properly chosen 
Boolean function fix), as explained in the previous sec- 
tions. 

Further, towards inverting the phase, we will use an- 
other symmetric Boolean function g{x), different from 
f(x), where g(x) = 1, when wt(x) = w, and g(x) = 0, 
otherwise. Based on g(x), we implement the inversion 
operator as O g , that inverts the phase of the states \x) 
where {x £ {0, 1}™ \wt(x) = w}. Thus, we consider the 
operator 

G t = [(2\y)(*\-I)O g f 

on 1*) to get |\l/ t ). 

Consider the n-qubit state \ty) = J2 s es Us \ s ) + 
Sse-tcil-xs -1 '*! 5 )' where u s ,v s are real and J2 s£S u 2 + 
Y^se{o,i} n \S v s = 1- For brevity, let us represent \^) = 
E u s \s) + £ v s \s) = u\X) + v\Y). That is, 

sGS sS { ,l}"\S 

u 2 = E u s an d y2 = E As- 
ses' se{o,i}«vs 

Theorem 4 Let u = sin (9, v = cos 6. The application 
of [(2|*)(*| - /)£><,]* operator on produces |* t ), in 
which the probability amplitude of \X) is sin(2i + 1)9. 

Proof: For t = 1, one can check that = [(2|*)(*| - 
T)Og]\V) = [(2|tf)(tf|)Cytf) -Cy*). Now substitut- 
ing the values of u,v, we get that = s'm39\X) + 
cos39\Y). 

Now we will use induction. Let the application of 
[(2|^)(*| - I)OgY operator on |*) updates the prob- 
ability amplitude of \X) as sin(2i6* + 9), for t = k. 
From the assumption we have [(2|*)(*| - I)O g ] k \^) = 
sm(9 + 2k9)\X) + cos(6> + 2k9)\Y). Now, for t = k + 1, 
it can be checked that [(2|*)(*| - i")0 s p +1 )|*) = 
sin(0 + 2(fc + 1)9)\X) + cos(9 + 2(k + 1)9)\Y). Thus, 
the proof. ■ 

We will now use such states \^t) in parity measure- 
ment. Consider that after the Deutsch-Jozsa algorithm 
we obtain a symmetric n-qubit state (before the mea- 
surement) 

\S n ) = ^2 a wt{x) \x), 

X£{0,1}*> 

such that (™)|au,| 2 = -7=, for some constant c. Thus, 
we have the initial amplitude of target states, {x € 
{0,l} n \wt(x) = w}, is sinfl = («)\a w \ = P^. For 
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n 4 


w — > 


l 


2 


3 


4 


5 


6 


7 


8 


4 


2>'/|0> 




0.833609 


0.981763 


0.833609 


- 


- 


- 


- 


- 




/ 




01 


02 


05 


- 


- 


- 


- 


- 




r 




0.468136 


0.298698 


0.468136 


- 


- 


- 


- 


- 






<)n 


0.5625 


0.375 


0.5625 














® 


0.421875 


0.375 


0.421875 












5 


v' f \o) 




0.748304 


0.92852 


0.92852 


0.748304 


- 


- 


- 


- 




f 




03 


02 


05 


16 


- 


- 


- 


- 




r 




1.42458 


0.313077 


0.313077 


3.57542 


- 


- 


- 


- 






§n 


0.703125 


0.625 


0.625 


0.703125 


- 


- 


- 


- 




® 


0.4096 


0.3456 


0.3456 


0.4096 










6 






0.730278 


0.823495 


0.954987 


0.823495 


0.730278 


- 


- 


- 




/ 




03 


02 


05 


OA 


29 


- 


- 


- 




r 




1.48129 


0.357282 


0.277975 


0.357282 


4.51871 


- 


- 


- 






§n 


0.585938 


0.527344 


0.3125 


0.527344 


0.585938 


- 


- 


- 




m 


0.401878 


0.329218 


0.3125 


0.329218 


0.401878 








7 


2?'/|0) 




0.704306 


0.754753 


0.907588 


0.907588 


0.754753 


0.704306 


- 


- 




/ 




07 


60 


05 


OA 


53 


4A 


- 


- 




r 




2.44507 


5.93733 


0.27984 


0.27984 


5.93733 


2.44507 


- 


- 




25/10) s 




0.683594 


0.512695 


0.546875 


0.546875 


0.512695 


0.683594 


- 


- 




m 


0.396569 


0.318745 


0.293755 


0.293755 


0.318745 


0.396569 






8 


V' f \Q) 


xjn 


0.698181 


0.710643 


0.813922 


0.92625 


0.813922 


0.710643 


0.698181 


- 




f 




3F 


CO 


BF 


AO 


AF 


AC 


AD 


- 




r 




5.51859 


6.91248 


7.69903 


7.74472 


7.69903 


6.91248 


5.51859 


- 






<)n 


0.598145 


0.553711 


0.413574 


0.273438 


0.413574 


0.553711 


0.598145 


- 




m 


0.392696 


0.311462 


0.281632 


0.273438 


0.281632 


0.311462 


0.392696 




9 


V' f \Q) 




0.684842 


0.651002 


0.76886 


0.884277 


0.884277 


0.76886 


0.651002 


0.684842 




f 




OF 


180 


OD 


140 


15F 


6A 


153 


16A 




r 




3.4566 


7.86171 


0.858163 


8.7469 


8.7469 


0.858153 


7.86171 


3.4566 








0.672913 


0.430664 


0.415283 


0.492188 


0.492188 


0.415283 


0.430664 


0.672913 




m 


0.389744 


0.306102 


0.273129 


0.260182 


0.260182 


0.273129 


0.306102 


0.389744 



TABLE I: Probability values using biased Hadamard transform (in this case we provide the corresponding symmetric function 
/ represented as a hexadecimal number of the (n + 1) length bit-string /„, f n -i, ■ ■ ■ , /i, fo and the value of r), using (standard) 
Hadamard transform and the method of [J]. 



large n, one can approximate it as 9 = -j^= and hence 
we need t iterations of Grover like strategy such that 

One important difference with Grover algorithm 
here is that we have good (almost exact) estimate 
of t, which is not known priori for application in 
search algorithms. After the application of Grover like 
strategy, we will get another symmetric n-qubit state 
\T n ) = E < t( Jx) such that OU 2 will be very 

z6{0,l}" 

close to 1 and the parity measurement will produce \D™} 
mostly in one step with very high probability. Thus 
the exact strategy is similar to Algorithm 1 (Algorithm 
2 can be modified with a similar way) in the previous 
section, where we add one more step as follows. 



Algorithm 3 

1 . Let / G SB n be as explained in Theorem [1] to maxi- 
mize the Walsh spectrum values at weight w. 

2. Use any of the above three strategies (our strategies 
exploiting Hadamard or biased Hadamard transform or 
the strategy of [l|) to obtain a symmetric n-qubit state 
\S n ) = E a wt(x) \x), such that C)\a w \ 2 is 

xG{0,l}" 

3. Use Gt on \S n ), t many times, where t is O(^fn) to 

obtain |T«) = E <t( x )\ x ) such that OKI 2 is 
ze{o,i}" 

very close to 1. 

4. Use parity measurement strategy. If the ancilla state 
is measured at the basis U W \Q, then \D™) is successfully 
obtained. Else go to Step 2 and iterate. 
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Example 3 Let us consider n = 8, w = § = 4. One may 
note that using our technique in Section \IIID\ or using 
the technique o/|3/, one can obtain \D™) with probability 

%^ = 0.2734375. Thus, sm9 = 0.523. Taking, 

t = 1. one may note that sin(2i + 1)9 = sin W ps 0.9968. 
Hence after the application of Grover algorithm with t = 
1, the probability of obtaining |-D") increases to 0.9968. 

B. Potential Advantage 

The advantage of this method over Algorithm 1 is as 
follows. First, we need O(-tfn) steps using Grover algo- 
rithm in each run and then a parity measurement should 
provide Thus, we get a quadratic speed-up (which 

is quite natural) over just using Dcutsch- Jozsa algorithm. 
Second, the number of parity measurement was 0(y/n) 
in the earlier case, once in each iteration. Here it is only 
a very few (may be 1 in most of the cases). 

It is natural to use Grover algorithm for amplitude 
amplification, but in the known applications, the num- 
ber of target states for which the amplitude is increased 
are not large. That guarantees the efficient implementa- 
tion of the phase reversal circuit. In the case of preparing 
l-D"), the situation is different as we need to amplify the 
amplitude at (™) many points and this could be expo- 
nential. Thus, it is an important question whether the 
phase reversal can be implemented efficiently. In this 
case, this can be achieved as the phase reversal can be 
implemented with symmetric functions, the implementa- 
tion of which is efficient [||. Further, as the amplitude af- 
ter the Deutsch- Jozsa algorithm with (biased) Hadamard 
transform is known in this case, the number of iterations 
in Grover algorithm will be accurately estimated apriori 
and thus the problem of obtaining the proper stopping 



criteria will not occur in this case. 



VI. CONCLUSION AND OPEN PROBLEMS 

In this work, we study several quantum algorithms to 
construct arbitrary Dicke state in a disciplined manner. 
The key idea is to find a suitable symmetric Boolean 
function for Dcutsch- Jozsa algorithm for the given n 
and w, use of the Grover algorithm and the general- 
ized parity measurement strategy. Further, we show 
that it is possible to obtain improved results using bi- 
ased Hadamard transform suitably. Our results improve 
the probabilities obtained in [l[ and thus provide faster 
method to construct Dicke states. The problem open 
in this area is to characterize the enumeration results 
in case of modifying the Deutsch- Jozsa algorithm with 
biased Hadamard transform. Obtaining the exact bias 
(^) in biased Hadamard transform with the correspond- 
ing symmetric function to optimize the probability cor- 
responding to the Dicke state seems to be an interesting 
problem. 

Though we look at the problem from theoretical an- 
gle, the algorithmic blocks used by us have experienced 
major advancement towards actual implementation. One 
may refer to flol Section 7] for literature related to im- 
plementation of quantum gates as well as Deutsch- Jozsa 
algorithm, Grover algorithm and several measurement 
techniques. As example, the idea of implementing bi- 
ased Hadamard transform is related to the Fabry-Perot 
cavity [H Page 299]. 
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